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1 INTRODUCTION 



Standard inflationary models (e.g. Peacock 1999 ) predict 
that the large-scale structure of the present day Universe 
originates from primordial density perturbations amplified 
by gravitational instability. Under the influence of its own 
gravity, the initial Gaussian random field of the density fluc- 
tuations evolves into a strongly non-Gaussian density field. 



ABSTRACT 

We analyse the redshift space topology and geometry of the nearby Universe by com- 
puting the Minkowski functionals of the Updated Zwicky Catalogue (UZC). The UZC 
contains the redshifts of almost 20,000 galaxies, is 96% complete to the limiting mag- 
nitude mz w = 15.5 and includes the Center for Astrophysics (CfA) Redshift Survey 
(CfA2). From the UZC we can extract volume limited samples reaching a depth of 
70/i _1 Mpc before sparse sampling dominates. We quantify the shape of the large- 
scale galaxy distribution by deriving measures of planarity and filamentarity from the 
Minkowski functionals. The nearby Universe shows a large degree of planarity and 
a small degree of filamentarity. This quantifies the sheet-like structure of the Great 
Wall which dominates the northern region (CfA2N) of the UZC. We compare these 
results with redshift space mock catalogues constructed from high resolution A-body 
simulations of two Cold Dark Matter models with either a decaying massive neutrino 
(rCDM) or a non-zero cosmological constant (ACDM). We use semi-analytic mod- 
elling to form and evolve galaxies in these dark matter-only simulations. We are thus 
able, for the first time, to compile redshift space mock catalogues which contain galax- 
ies, along with their observable properties, rather than dark matter particles alone. 
In both models the large scale galaxy distribution is less coherent than the observed 
distribution, especially with regard to the large degree of planarity of the real survey. 
However, given the small volume of the region studied, this disagreement can still be 
a result of cosmic variance, as shown by the agreement between the ACDM model and 
the southern region of CfA2. 

Key words: methods: numerical; methods: statistical; galaxies: statistics; large-scale 
structure of Universe 



tities only describe a Gaussian density field in a unique and 
complete way. 

One possible approach to capturing higher-order fea- 
tures of the galaxy distribution involves geometrical descrip- 
tors. A number of these statistics have been suggested to 
date. Among others, recent morphological analyses of the 
nearby Universe have been performed through measu ring 



Therefore, in order Lo describe the observed galaxy disXribu- 
tion qu antitatively, we need statistics beyond the traditional 
Fourier transform pair P(k) and £(r), the power spectrum 
and the two-point correlation function. In fact, these quan- 

3 Present address: Universita degli Studi di Torino, Diparti- 
mento di Fisica Generale "Amedeo Avogadro" , Via Pietro Giuria 
1, 1-10125 Torino, Italy. 



the void probabili ty function of the CfA2 catalogue (Vo- 



geley et al. 1994a), ap plying percolatio n techniques to the 



IRAS 1.2Jy catalogue (Yess et al. 1997) and the Las Cam 



panas Redshift Survey (LCRS) (^handarin fc Yess 1998), or 



through various types of shape statistics, as the moment 



of-inertia meth od suggested by Babul fc Starkman (1992) 
and applied by Dave et al. (1996)|to the CfAl redshift sur- 



vey and by sathyaprakash et al. (1998) to the IRAS 1.2Jy 
catalogue. The most wide-spread technique is probably 
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Figure 1. Simple structures with positive (eg 
characteristic V3. 



left panel), vanishing (loop, middle panel), and negative (double loop, right panel) Euler 



the genus statistics introduced by Gott III et al. (1986)1 



which has been used on various IRAS catalogues 
et al 



1)92 



Protogcros fc Weinberg 1997; 



Canavezes et al 



199cj[f Jringel et al. 1998Q , on the CfA2 redshift survey dVo- 



geometric quantity 



V u 



geley e 1 al. 19 94b[) an ~ in a modifi ed way, on the nearly 
two-dimensional LCRS (dolley 1997). With galaxy clusters, 



scales of several hundred h Mpc have been probed using 
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different morphological statistics (Coles et al. 1998; 



et al. 1 )92 ) , including the Minkowski functionals ( Kerschei 
et al. 1)97). 



Plionis 



Table 1. Minkowski functionals V M 
corresponding geometric quantities. 



expressed in terms of the 



Among this plethora of methods, Minkowski function- 
als currently represent the most comprehensive description 
of the topological and geometric properties of a three- 
dimensional structure. Minkowski functionals include both 
the percolation analysis and the genus curve. Minkowski 



spectively wit hout discontinuities, we c an address this point 



here. In fact, 



Geller fc Huchra (1989) 



were the first to ad- 



functio nals can be a pplied both to density fiel ds (Schmalz- 
ing & I uchert 1997) and to discrete point sets (Mecke et al 



1994), where the latter approach has the advantage of intro- 
ducing just one single diagnostic parameter (see below for 
details), rather than, for example, a smoothing length and a 
threshold as in the genus statistics. Because Minkowski func- 
tionals are additive they are robust measures even for sparse 



dress the two-dimensionality of the structures by identifying 
the "Great Wall" in the CfA2 redshift survey. On the other 
hand, with the UZC, we can only probe the topology of a 
volume 70/i~ 1 Mpc deep, which is too small to be represen- 
tative of the Universe. 

We compare the UZC redshift space topology with mock 
galaxy redshift surveys extracted from A-body simulations 
of a Cold Dark Matter (CDM) model i ncluding galaxies 



modelled with semi-analytical techniques ( Kauffmann et al 



cosmological models efficiently ( schmalzing et al. 1999a ) 
Finally, Minkowski functionals can be used to define a pla- 
narity V and a filamentarity T and t hus quantify our intu- 
ition about the shape of a structure (Sahni et al. 1995). 

Here, we study the redshift space topology of the nearby 
Universe by analys ing the Updated Zwicky Catalogue (UZC, 
Falco et al. 1999) which contains the redshifts of almost 



sampletj. Minkowski functionals also_ discriminate dittcrcnt 1999^. Therefore the density bias, the distribution of galax 



20,000 galaxies and is 96% complete in more than one fifth 
of the sky to the limiting Zwicky magnitude m z w = 15.5. 
The UZC includes the CfA2 Redshift Survey ([Huchra ct 
al. 199 fl; |Huchra et al. 1995|; [Huchra et al. 199C|; |Gellei 



Hucjira 1989|; [Hu chra et al. 1983; de Lapparent et al 
1982|). In our analysis, we use the 



1986 ; 1 )avis et al. 



Minkowski functionals and their associated shape quanti- 
fiers planarity V and filamentarity T which we define in 
Appendix [Bj. Recently, a similar analysis on the nearby Uni- 



ies relative to the dark matter, is taken self-consistently into 
account by including the physics of galaxy formation and 
evolution explicitly. 

In Section Q we review our analysis method. Section H 
briefly describes the properties of the UZC catalogue and 
presents the results of our topological analysis. We com- 
pare the northern region of the CfA2 redshift survey with 
the mock catalogues in Section r&. We summarise in Sec- 
tion o. Two appendices provide important but gory details 
of Minkowski functional analysis. 



2 METHOD 



2.1 Minkowski functionals 



verse h as b een performed by using the LCRS (Bharadwaj 
et al. 1999). The main advantage of using the UZC over 



the LCRS is the geometry of the survey: the LCRS consists 
of six separated slices only 1?5 thick; thus, an analysis of 
the full three-dimensional structu re cannot be performed. 
In fact, Bharadwaj et al. (1999) are unable to discrimi- 



nate whether the "filaments" seen in the LCRS are actu- 
ally sections of "sheets" . Because the UZC consists of two 
coherent regions covering ~ 14% and ~ 7% of the sky re- 



Minkowski functionals are named after Minkowski (1903) 
in acknowledgment of his pioneering contributions to inte- 
gral geometry (fo r an up-to-dat e overview of the the ory 
see, for example, Schneider 1993). Mecke et al. (1994) in- 
troduced the Minkowski functionals into cosmology as de- 
scriptors of the geometry and topology of a distribution of 
points. One appealing property of the Minkowski functionals 
is their uniqueness and completeness as morphological m ea- 
sures, which follows by the theorem of Hadwiger (1957) on 
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Figure 2. Surface area and the integrated mean curvature, 
namely the Minkowski functionals Vi and V2, for spheroids of 
varying axis ratio. Prolate spheroids ("cigars") are characterised 
by large surface area and comparatively small mean curvature, 
while the situation is reversed for oblate spheroids ("pancakes"). 



very general requirements. Furthermore, the four Minkowski 
functionals in three dimensions can be interpreted as well- 
known geometric quantities, as summarised in Table DJ. 

Figure |lj shows examples of bodies with different Euler 
characteristic V3. The Euler characteristic is a purely topo- 
logical quantity and measures the connectivity of a set. It 
is equal to the number of parts minus the number of holes. 
Hence the egg (left panel) has V3 — 1, while the single (mid- 
dle panel) and double loop (right panel) have V3 = and 
Vz — — 1, respectively. The behaviour of the surface area Vi 
and mean curvature V2 is illustrated in Figure ti on a simple 
family of sph eroids, whose M inkowski functionals are known 
analytically (Hadwiger 1955). By varying the axis ratio of 
the spheroids while keeping the volume fixed, we obtain a 
curve in the V1—V2— plane. Values are given in units of the 
corresponding functionals of a sphere of the same volume, so 
we are indeed showing the effects of shape. A sphere would 
lie at the cusp at (1,1). Prolate spheroids ("cigars") are 
characterised by large surface area and comparatively small 
mean curvature, while the situation is reversed for oblate 
spheroids ( "pancakes" ) . 

We consider a volume-limited subsample from a galaxy 
catalogue as a set of points {x;}, where i — 1 . . . N. Since the 
geometry of single points is trivial, we decorate each point 
with a ball B of radius r, and measure the Minkowski func- 
tionals of the resulting union set, the so-called Boolean grain 
model (Wicksell 1925). In short, we construct the union set 



A r = (J B Xj 



(1) 



Figure 3. The volume fraction of the reduced window as a func- 
tion of the shrinking distance. The dashed line corresponds to full 
sky coverage of 4tt; the solid line is for the CfA2n window. 



radius r. By increasing the size of the balls, connections are 
established between points farther apart, until the survey 
volume is completely filled by the Boolean grain model. The 
radius can therefore be used as a diagnostic parameter. 

2.2 Incomplete sky coverage and finite survey 
depth 

Currently, estimation of geometrical quantities such as 
Minkowski functionals relies on homogeneous sampling of 
the underlying point distribution. Galaxy catalogues, how- 
ever, are usually magnitude-limited, and therefore intrinsi- 
cally inhomogeneous. The standard solution is to use a series 
of volume-limited subsamples. However, even with volume- 
limited subsamples, care needs to be taken to obtain mean- 
ingful results that are not affected by the sample geometry. 
We summarise our approach here. An extensive description 
can be found in Appendix L\L 

For calculating the volume functional wo, we use the 
unbiased estimator 



Vo 



V (A r n D r ) 
V (D r ) ' 



(2) 



where D r denotes the part of the volume, or window, D that 
is farther than r from the boundary. For all other Minkowski 
functionals, we use minus estimators constructed from par- 
tial Minkowski functionals. An unbiased estimator of the 
volume densities w M of the Minkowski functionals is obtained 
through 



X^r> 2r (xi)^(A r ; Xi ) (jt= 1,2,3), (3) 



\D 2r 



and look upon the quantities V^ (A r ) as functions of the ball 
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where Z?2r denotes the reduced sample window, and XD 2r IS 
its characteristic function. 

This approach works for arbitrary sky coverage, includ- 
ing the non-convex boundaries of the CfA2 survey. More- 
over, at any radius, we use all available points, since points 
outside the reduced sample window still contribute as neigh- 
bours of points remaining inside. The only d rawback of our 
approach, when compared to the method of |3chmalzing et 
al. (1996)1 is that the window now shrinks by twice the ra- 
dius. This fact considerably reduces the maximum allowed 
size of the Boolean ball. 

We may quantify the loss of accuracy by the ratio 
|-D r |/|-D| between the volumes of the reduced and the origi- 
nal windows. For complete sky coverage, we have 



\JM 

\D\ 



(>-i) 



= 1- 



(4) 



where R gives the depth of the volume-limited sample and r 
is the distance to be kept from the boundary. If the sky cov- 
erage is not complete, this fraction decreases rapidly. This 
is illustrated in Figure H, which compares the reduced win- 
dow size for full sky coverage, and for the 10% sky coverage 
of the CfA2n window (defined in Section |f.l| ) . Apparently, 
analysing the CfA2 catalogue with Minkowski functionals 
becomes impractical at radii above 10/i -1 Mpc. We set our 
maximum Boolean ball radius at 5/i _1 Mpc, so our results 
are not affected by this problem. 



3 THE UPDATED ZWICKY CATALOGUE 

We now briefly describe the pro pert ies of the Updated 
Zwicky Catalogue (UZC, Section |3.l|) . In Section 3.2 we 
compute its Minkowski functionals with the method based 
on the Boolean grain model. 



3.1 Description 

Falco et al. (1999) provide new measurements of the posi- 
tions and redshifts of most galaxies in the Zwicky Catalogue, 
which has been the basis for the CfAl and CfA2 redshift sur- 



veys (see Falco et al. 199£ and references therein) . The UZC 
represents a homogeneous and well calibrated redshift cata- 
logue of almost 20,000 galaxies. Furthermore, the catalogue 
is now publicly available. 

The UZC covers the northern celestial hemisphere 
which is not obscured by the Galactic disk; the redshift cat- 
alogue is ~ 96% complete to a limiting Zwicky magnitude 
mzw = 15.5. 

The CfA2 redshift surveys have been compiled over 
many years and this name has indicated progressively larger 
areas of the sky, as the surveys were reaching completion. 
The UZC contains all of the CfA 2 survey analysed i n the lit- 
erature to April 1999. Following Falco et al. (1999), we refer 



to the area with right ascension ranges 8 < 01950 < 17 11 and 
20 < «i95o < 4 h and declination range — 2? 5 < £1950 < 50° 
as the CfA2 region. The first right ascension range defines 
the CfA2N region containing 12,082 galaxies with measured 
redshift. The regions between 3 h and 4 h and between 20 h and 
22 h are sparsely populated in the catalogue largely because 
of galactic obscuration. Therefore we restrict the CfA2S re- 
gion to the 4,436 galaxy redshifts in the right ascension 




Figure 4. Cone diagram of the volume— limited subsample of 
70h~ Mpc depth we analyse here from the UZC. The larger and 
smaller cones correspond to the northern and southern region, 
respective ly. The large do t marks the location of our galaxy and 



alludes to Adams (1986) 



range 22 1 < 01950 < 3 . Both areas used in our analysis are 
more than 98% complete. Moreover the UZC contains some 
galaxies fainter than the magnitude limit. These galaxies 
are mainly in multiplets which were unresolved by Zwicky. 
They affect the properties of the galaxy distribution on very 
small scales but should leave our analysis unchanged. Note 
that the actual number of available redshifts in the UZC is 
currently increasing, as the catalogue is being updated. 

Visual inspection shows well-defined structure in the 
CfA2 surveys, with large voids and sheet-like structures. 
Figure W shows a cone diagram of the volume-limited sam- 
ples from CfA2N and CfA2S we analyse here. The dense 
feature close to the boundary of the northern part of the 
catalogue is the "Great Wall" (feller fc Huchra 1989, see 



also Figure 6 of Falco et al. 1999). In CfA2S, the apparent 
galaxy concentration is the Perseus-Pisces supercluster. 

3.2 Topology and geometry 

In our analysis, we use volume-limited samples 70/i _1 Mpc 
deep. This provides a good compromise between object num- 
ber (of all volume-limited subsamples, the 70/i _1 Mpc one 
includes the highest number of objects) and depth (parts 
of the Great Wall are closer than 70/i _1 Mpc). Since the 
scaling behaviour of Minkowski functionals with the num- 
ber density is not known in general, we take care to always 
analyse samples with the same number density of galaxies. 
From all volume-limited samples, we choose a set of slightly 
smaller random samples such that the number density cor- 
respond s to exactly 1,000 galaxies in the CfA2n volume (see 
Section 4.1 for a definition). 
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Figure 5. Minkowski functionals of CfA2N (solid line) and CfA2S (dashed line). All curves were calculated from volume— limited samples 
70fe -1 Mpc deep. No te that the CfA2N Minkowski functionals are systematically larger than the CfA2S ones, indicating larger clustering 
(Marzke et al. 1995h. 



The results of our analysis of the UZC are shown in Fig- 
ure pi. Instead of the Minkowski functionals w M themselves, 
we plot the quantities 1 — V^/V^B). A thorough motiva- 
tion of their use is given in Appendix |B|; most notably they 
are exactly equal to zero for a stationary Poisson process, 
i.e. a random distribution of points with spatially constant 



density. Generally speaking, values above zero reveal a clus- 
tered distribution of galaxies. Note, however, that the V M 
can be written as an alternati ng series containing the h ier- 
archy of correlation functions ( jchmalzing et al. 1999b ), so 
clustering does not simply mean that the two-point func- 
tion is positive. As an example, consider 1 — Vo/Vo(B). By 



© 1999 RAS, MNRAS 666, 1-14 



6 Jens Schmalzing and Antonaldo Diaferio 



Equation (pl|) this quantity increases monotonica lly with 
1 — Vo, the void probability function (White 1979). Hence 



larger values in our plots indicate larger voids in the point 
distribution, which in turn means that points have to cluster 
more tightly. 

Some of the functionals, most notably the volume vo, in- 
dicate stronger clustering for CfA2N than for CfA2S. In fact, 
CfA2N contains more clusters than CfA2S and the Great 
Wall is more apparent than the Perseus-Pisces superclus- 
ter in CfA2S. The two-point correlation function already 
showed that the CfA2N clu stering is peculiarly large when 



Diaferio et 



compared to other surveys (Marzke et al. 1995 
al. 1999^). 

We use a toy model to characterise the geometry of 
CfA2: we assume that galaxies are distributed in one- 
dimensional filaments, two-dimensional sheets or in a homo- 
geneous field. We can compute the Minkowski functionals for 
this model analytically. Two coefficients, T and V, give the 
relative contribution of filaments and sheets in the galaxy 
distribution (see Appendix M for details). Note, however, 
that this approach has the main shortcoming of assigning 
structures which are loosely clustered to the planarity com- 
ponent, because alternative assignments would yield consid- 
erably worse fits. Therefore, T and V yield the relative im- 
portance of the two structural components, but their exact 
values should be considered cautiously. 

By fitting our toy model to the CfA2N Minkowski func- 
tional curves, we obtain T — 0.12 and V = 0.67 for the 
filamentarity and planarity contribution, respectively. The 
high degree of planarity is expected from the visual appear- 
ance of the Great Wall in the catalogue. For the CfA2S we 
find T = 0.01 and V = 0.75. The large value of V in both 
regions reflects the tendency of galaxies to be distributed 
in two-dimensional structures or in loose structures rather 
than in filaments. Moreover, the smaller size of the CfA2S 
region also tends to increase slightly the signal for planar 
and loose structures. 

Because of its simplicity, our toy model provides a 
straightforward way of measuring planarity and filamentar- 
ity, but it is inadequate to model more complex realistic 
situations. We are planning some refinements in order to 
obtain results in better agreement with our intuitive inter- 
pretations of the galaxy distribution. 



4 MOCK GALAXY REDSHIFT SURVEYS 




Figure 6. Cone diagrams of two mock catalogues. The upper and 
lower cones are a volume-limited subsample constructed from the 
CfALF ACDM model and the CfALF -rCDM model, respectively. 
The observer galaxy is indicated by the large dot. 



Table 2. GIF Simulations 



model 


tto 


n A 


h 


°"8 


nip 


L 


ACDM 


0.3 


0.7 


0.7 


0.90 


1.1 


111 


rCDM 


1.0 


0.0 


0.5 


0.60 


1.0 


85 



Parameters of the two GIF simulations ( Kauffmann ct al. 1999a ) 
used in this paper. The Hubble constant h, the particle mass 
m p , and the comoving size L of the simulation box are in units of 
H = 100kms~ 1 Mpc _1 , 10 10 /i _1 Mq, and h~ 1 Mpc, respectively. 

a% of the power spectrum of the initial density perturbations 
was chosen to reproduce the present day abundance of rich 
galaxy clusters. Table H summarises the parameters of the 
models. 



The s imulations were run with Hydra ( Pearce fc Couch 



man 1997J ) . the parallel version o f the AP M code (Couch- 
1991|; ICouchman et al. 199q) , kindly provided by Jenk 



We co mpare the morphology of the CfA2n r edah ift survey 



ins et al. 



(1997)1 



with two variants of a CDM model. In Section 4.1 we briefly 



describe the simulatio ns an d how we model galaxy formation 
and evolution. Section 4.2 discusses how the luminosity limit 



of the galaxy sample, the shape of the sample volume and 
the peculiar velocity distortion s in redshift space affect the 
Minkowski functionals. Section 4.3 contains the comparison 



of CfA2n with the mock galaxy redshift surveys. 



4.1 Models 



We consider two sim ulations from the GIF project ( Kauff- 



mann et al. 1999a ): a CDM model with either Oo = 1 
(rCDM) or Cl = 0.3 and A = 0.7 (ACDM). The simulations 
contain 256 3 dark matter particles each. The normalisation 



rvauttmann et al. (tyyyaj combine these high resolution 
A^-body simulations with semi-analytic modelling to form 
and evolve galaxies within dark matter haloes: gas cooling, 
star formation, supernova feedback, stellar evolution and 
galaxy-galaxy merging are the relevant processes included. 
At any cosmological epoch, the observable properties of the 
galaxies, namely luminosity, colour, star formation rate, and 
stellar mass, are the result of the galaxy's merging history, 
its interaction with the environment and the passive evolu- 
tion of its stellar content. 

The models predict photom etric and clustering pro per 
ties of galaxies at re dshift z = ( Kauffmann et al. 1999a ; Di 
aferio et al. 1999a), and their evolution to hi gh redshift 



( [Kauffmann et al. 1999b ; Diaferio et al. 1999b) which are 
in reasonable agreement with observations. Moreover, these 
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analyses suggest further diagnostics for constraining the cos- 
mological model and determining the relevance of the dif- 
ferent galaxy formation processes. 

These simulations represent the only models, available 
to date, where galaxies form and evolve in a self-consistent 
cosmological model and within a volume of the Universe suf- 
ficiently large for an analysis of the large scale galaxy distri- 



bution (see also Pearce et al. 1999). In fact, Diaferio et al 



(1999a) extract mock galaxy redshift surveys from the sim- 
ulations at z — for a direct comparison with the CfA2N 
redshift survey, whose volume ~ 7 x 10 5 /i _3 Mpc 3 is com- 
parable to the volume of the simulation boxes. Specifically, 
they consider the area 8 h < a 1950 < 17 h , 8?5 < 5i 950 < 44?5 
which has a declination range smaller than CfA2N. We refer 
to this region as the CfA2n region. In order to resemble the 
CfA2n appearance in the mock catalogues as closely as pos- 
sible, the mock redshift surveys were compiled by locating 
the observer home galaxy within the simulation box on a 
galaxy similar to the Milky Way at ~ 70/i _1 Mpc away from 
a Coma-like massive cluste r. 

Diaferio et al. (1999a) analyse both the CfA2n and the 
mock redshift surveys with the same standard observational 
techniques. They show that the small scale clustering prop- 
erties of galaxies, namely galaxy group properties, redshift 
space correlation function and pairwise velocities, are repro- 
duced by both models independently of fio- Despite this 
success, a visual inspection of the mock redshift surveys 
shows remarkable differences in the large scale distribution 
of galaxies. Both models fail to produce structures as coher- 
ent and as sharply defined as in the CfA2n survey, although 
ACDM catalogues seem to be in better qualitative agree- 
ment with the data (see their Figures 5-8). Two examples 
of volume-limited subsamples constructed from these mod- 
els are shown in Figure q. 

The B-band galaxy luminosity function predicted by 
the semi-analytical models is not a good fit to the CfA 
survey luminosity function ( Marzke et al. 1994 ). To inves- 
tigate the effect of the lum inosity function on the small- 
scale clustering properties, Diaferio et al. (1999a) assign 
new luminosities to the galaxies in order to reproduce the 
CfA luminosity function exactly, while preserving the galax- 
ies' luminosity rank. We thus have two sets of mock cata- 
logues, where the galaxies have luminosities according to 
the semi-analytical model luminosity function (SALF) or 
the CfA luminosity function (CfALF). We have analysed 
both the SALF and CfALF mock catalogues. Similarly to 
other measures related to the underlying mass distribution, 
the Minkowski functionals are basically independent of the 
luminosity function adopted. Throughout this analysis, we 
show results from the CfALF catalogues alone. Moreover, we 
consider ten mock catalogues for each simulation to assess 
the robustness of our results. 



4.2 Luminosity, geometry and redshift space effect 

Kauffmann et al. (1999a) show that the two-point corre- 
lation function is similar for galaxy samples with different 
magnitude limits. The Minkowski functionals are able to ex- 
tract more information from the galaxy distribution. In fact, 
bright galaxies cluster stronger than faint ones. Figure Q 
shows the Minkowski functionals for the ACDM simulation; 
rCDM yields very similar plots. In order to keep constant 



the statistical significance of different galaxy samples, we 
compute the Minkowski functionals only of subsamples con- 
taining the same number density of galaxies, and average 
over ten subsamples. The solid lines show the Minkowski 
functionals for galaxies brighter than Mb = —17.5 + 5 log h. 
The Minkowski functionals for galaxies brighter than Mb = 
— 18.7 + 51og/i are indicated by dashed lines. This ab- 
solute magnitude corresponds to an apparent magnitude 
tub — 15.5 at 70ft -1 Mpc. The difference between the two 
curves is significantly larger than the scatter between differ- 
ent subsamples from the same set of galaxies. 

We have also checked that the Minkowski functionals 



computed with our procedure (Section 2.2) are insensitive 
to the shape of the sample volume. We do this by calculating 
the Minkowski functionals of the real space position of the 
galaxies averaged over the volume-limited subsamples ex- 
tracted from the ten mock redshift surveys. The dashed lines 
in Figure M lie within the scatter of the resulting Minkowski 
functionals (shown as the dark shaded area in Figure pf). 
Since they come from volumes which have differing shape 
(periodic, rectangular box versus CfA2n volume ) but con- 
tain galaxies with roughly the same luminosity, we can de- 
duce that the difference due to the sample geometry is negli- 
gible in comparison to the scatter between individual mock 
catalogues. Note that for mock catalogues, the scatter is 
much larger than for subsamples from the whole simulation 
box, simply because a mock catalogue contains fewer galax- 
ies. 

Redshift space distortions tend to destroy structure on 
small scales, because groups and clusters in real space are 
diluted in redshift space by the finger-of-god effect. Figure w 
shows this difference in the ACDM model. The Minkowski 
functionals of the real and redshift space galaxy distribu- 
tions almost coincide at larger values of the Boolean ball 
radius, because redshift space distortions become less rele- 
vant on large scale. 

4.3 Comparison with the CfA2n 

Figure |9| compares the morphological properties of CfA2n 
with the mock catalogues: neither of the two models is ca- 
pable of reproducing the CfA2n properties. At Boolean ball 
radii ~ 1 — 2/t -1 Mpc, all the four Minkowski functionals of 
CfA2n are five times the r.m.s. above the average mock cat- 
alogue of the ACDM. Given the rather small number of ten 
realisations, we can conclude that the CfA2n clustering is 
not reproduced by the ACDM model at the 90% confidence 
level. The rCDM model shows an even smaller degree of 
clustering. At larger radii and for the mean curvature V2 and 
Euler characteristic V3, the disagreement between the mod- 
els and observations is less dramatic, but also less significant, 
because the sample volume is filled by a considerably smaller 
number of Boolean balls. It is crucial to recall, however, that 
the scatter of quantities derived from the variance over all 
ten mock catalogues underestimates the sampling variance, 
because the mock surveys are constructed from the same 
small parent simulation. 

The failure of the ACDM in reproducing the CfA2n 
Minkowski functionals is particularly interesting, because, 
on average, this model fits other observations reas onably 



well and in any c a se better than the rCD M model (Kauff- 
mann et al. 1999^ ; JKauffmann et al. 1999b] ) . The Minkowski 
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Figure 7. Effect of luminosity and geometry on the Minkowski functionals of the CfALF ACDM model. The dashed and solid lines 
are for all galaxies in the simulation box brighter than M-g, = —18.7 + 51og?i and Mg = —17.5 + 51og/i respectively. The shaded area 
indicates the lc-spread over ten subsamples. 



functionals confirm that ACDM should be preferred to 
rCDM. 

It is worth noting that, because of the smaller degree 
of clustering of CfA2S when compared to CfA2N (see Fig- 
ure £3|), the difference between the CfA2S and the models is 
less pronounced, although still significant. The result is an 



exampl e of the power of the M inkowski functionals analysis. 
In fact, Diaferio et al. (1999a) show that the two point corre- 



lation function of the mock catalogues agrees quite well with 
the CfA2S function. As expected, the Minkowski function- 
als are able to extract more information on the clustering 
properties of galaxies. 
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Figure 8. Minkowski functionals of the ACDM mock catalogues in real space (dark shaded area) and in redshift space (light shaded 
area). The shaded area indicates the r.m.s. spread over the ten mock catalogues. 



We can also perform the fitting procedure described 
in Appendix |b| to estimate the planarity and filamentarity 
of the structure contained in the mock catalogues. ACDM 
yields T = 0.01 and V = 0.64, while rCDM results in 
T — 0.02 and V = 0.60 for the average mock catalogues. As 
pointed out in Section 3.2, the large values of V compared 
to T indicate that the galaxies tend to be distributed in 



two-dimensional or loose structures rather than filaments. 
Compared to the values determined from the real CfA2N 
catalogue in Section |J (T = 0.12 and V = 0.67), V is 
~ 10% smaller in the models than in CfA2N, again indicat- 
ing a smaller degree of clustering. 

We cannot exclude that the failure of the models in re- 
producing the Minkowski functionals of CfA2N is due to the 
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Figure 9. Structure in the CfA2n redshift survey, and in the mock catalogues extracted from the ACDM (dark shaded area) and the 
rCDM (light shaded area) models. The areas indicate the 1-cr scatter over the ten mock catalogues constructed for each model. All 
curves are for volume-limited samples 70h~ 1 Mpc deep. 



small size of both the simulation box and the real survey. 
In fact, the initial power spectrum of the density perturba- 
tions does not contain Fourier components with wavelength 
larger than the box size. This prevents the development of 
large-scale non-Gaussian features such as the Great Wall. 
Moreover, as already indicated by an analysis of the IRAS 



1.2Jy catalogue (Kerscher et al. 1998), a volume-limited 
subsample 70/i _1 Mpc deep is too small to be a fair sam- 
ple of the Universe. Therefore, our analysis is likely to be 
dominated by local effects, which should average out when 
deeper catalogues are considered. This hypothesis could al- 
ready be tested either with simulations of larger volumes 
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(see e.g. Doroshkevich et al. 1999), where one could esti- 



mate the cosmic variance, or with a constrained realisation 



(Mathis et al. 1999), which should not suffer from the fair 
sample problem. 



tions, but still the agreement between ACDM and CfA2 does 
not improve. 

Note that our Boolean grain model approach follows 
a direction complementary to the genus statistics an alysis 
generally used. The standard method introduced by Gott 



III et al. (1986) employs a density field smoothed on vary- 



5 CONCLUSION 



ing scales. Increasing the smoothing kernel width succes- 
sively er ases structure on small scales (ter Haar Romeny et 



de Lapjpenret ai. (WHO) ana juener | nucnra (iysy)\ nrst a l, 199l|). The analysis of |Vogeley et al. (1994b) 



for exam- 



pornted 1 ourcliat the ^ large-SCale galaxy dteWlbUtlOll ill the ple> employs widths between 6r'Mpc and 20/i Ll Mpc and 

addresses the issue of the Gaussianity of the distribution 
of the initial density pert urbations. We adopt the Boolean 
grain model approach by Mecke et al. (1994) and decorate 



CfA2 redshift survey has a two-dimensional sheet-like struc- 
ture rather than a one-dimensional filamentary appearance. 
Here we quantify this visual impression by analysing the red- 
shift space topology and geometry of the UZC. By construct- 
ing a Boolean grain model of balls centred around the galax- 
ies, we obtain unbiased estimates for the Minkowski func- 
tionals of the galaxy distribution. We then fit the Minkowski 
functionals to a toy model to extract quantitative informa- 
tion on the degree of filamentarity and planarity in CfA2. 
Wc find that filament- and sheet-like structures contribute 
12% and 67% to the Minkowski functionals of CfA2N re- 
spectively, while the remainder comes from a homogeneous 
distribution of galaxies. Our results therefore st rongly sug- 



gest that t he filaments identified in the LCRS (Bharadwaj 
et al. 1999) are actually sections of sheets. 



Apparently, our simple model is inadequate in other re- 
alistic situations, such as CfA2S (the southern part of the 
UZC), where looser structures are present. However, the en- 
couraging result we obtain with CfA2N indicates that it is 
worth pursuing this approach with some refinements. When 
deeper surveys as the Sloan Digital Sky Survey will be avail- 
able, our approach can represent a valuable tool to quantify 
the top ology of the galaxy di stribution. Most notably, con- 



trary to 5ahni et al. (1998) or Schmalzing ct al. (1999a) , the 
method presented here does not rely on shape measurements 
of isolated isodensity contours. Hence we can also obtain 
meaningful results in low to intermediate density environ- 
ments, where objects in a smoothed field already percolate 
and are inaccessible to measurements. 

We also compare CfA2n, a smaller volume of CfA2N, 
with mock redshift surveys extracted from a ACDM and a 
rCDM model. The models do not show the high degree of 
clustering of CfA2n. However, CfA2n seems to be quite a pe- 
culiar region of the Universe and the disagreement between 
the models and, for example, CfA2S is less dramatic. More- 
over, the Minkowski functionals of the ACDM are closer to 
observations than the Minkowski functionals of rCDM. In 
conclusion, although our topological analysis confirms ear- 
lier results which indicate that overall ACDM is a better 
representation of the real Universe than rCDM, we show, at 
the same time, that ACDM is not yet a satisfactory model 
of the Universe. 

Vogeley et al. (1994b) compare the genus statistics of 
CfA2 with a ACDM model. They find that this model is 
consistent with CfA2 on scales > 10fo -1 Mpc, but it fails to 
reproduce CfA2 on smaller scales. They use a simple bias 
scheme to identify galaxies and extract mock catalogues 
from A-body simulations. So the physics acting on small 
scales, which is missing from their simulations, might be re- 
sponsible for the disagreement. Here, we used semi-analytic 
modelling to form and evolve galaxies in A-body simula- 



each galaxy with a sphere of radius r varying from zero to 
5h~ Mpc. This procedure captures the morphology of dis- 
tinct structures on very large scales while preserving the 
small scale structure information. 

Finally, the small volume of both the survey and our 
simulation boxes represents a severe limit to our analysis, 
because we are unable to quantify the effect of cosmic vari- 
ance. In a forthcoming paper, we will use simulation boxes 
of 320/i -1 Mpc and 480/i _1 Mpc on a side to quantify these 
effects. We will also extrac t mock reds hift surveys of the 
Sloa n Digital Sky Survey (Gunn 1995) and the 2dF sur- 
vey (Maddox 1998). We will use a s imple biasing scheme to 



locate galaxies (e.g. Cole et al. 1996), rather than the semi- 
analytic modelling we have used here. With these larger vol- 
umes, we should be able to use our topological approach for 
discriminating between cosmological models successfully. 
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APPENDIX A: PARTIAL MINKOWSKI 
FUNCTIONALS 

Throughout our article, we consider the Minkowski func- 
tionals of a Boolean grain model A r constructed by placing 
balls of radius r around a set {x^} of N points contained in 
a spatial region D; 



\jB r 



(Al) 



For comparing samples of different size, we may define vol- 
ume densities u M of the Minkowski functionals of this objects 
by 

1 



\D 



V^Ar 



(A2) 
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if D is a large box with periodic boundary conditions, and 
\D\ denotes its volume. 

Unfortunately, for many practical applications, the sam- 
ple volume is not a box and does not feature periodic bound- 
aries. Therefore, one must use unbiased estimators for calcu- 
lating the Minkowski functionals; such a method of remov- 
ing the boundary effects given arbitrary sample geometries 
is provided through partial Minkowski functionals. 

First of all, to take into account the boundary effects 
on the volume functional vo, the sample window D must be 
reduced by the radius r to obtain an unbiased estimate. If 
D r denotes the set of all points that are farther than r from 
the boundary of the sample window D, we have 



Equation (A5) is already useful for the practical evalu- 



Vo 



V {A r r\D r 
V (D r ) 



(A3) 



Maurogordato & Lachieze-Rey (1987) introduced this esti- 



mator for calculating the void probability f unction of the 
CfAl catalogue. An extension was made by |5chmalzing et 
al. (199p)| to include all Minkowski functionals of a point set 
in a convex window. 

While for a smooth body, all Minkowski functionals v,x, 
p = 1,2,3 can be written as surface integrals, this is no 
longer possible for the Boolean grain model, which has edges 
and corners at intersections of two or more balls. Neverthe- 
less the contributions of these non-regular surface to the 
Minkowski functionals may still be evaluated using a more 
general con cept of surface integration th at includes edges 
and cor ners ( Federer 1959; Schneider 1993] ). Following Mecks 
et al. ( L994) , we arrive at the so-called partition formula, 
which sorts contributions to the Minkowski functionals from 
the various types of boundary, i.e. surface points belonging 
to one, two or three ballsQ The partition formula states that 

V^Ar) = 

N N N 

£ V^(A r ) + I £ V^\A r ) + \ £ V^ k \A r ). 



2 4V 6 

s,j=l i,j,k=l 



(A4) 



V M (A r ) denotes the contribution from the uncovered sur- 
face of the ball around point Xj, while V^ J {A r ) comes from 
the intersection line of balls i and j, and Vp i (A r ) is lo- 
cated at the corner made by the intersection of balls i, j and 
k. Rearranging the summations, we have 



VpiAr) = '£Vr(ArlXi) 



(A5) 



where the partial Minkowski functional V^(A r ;xi) sums up 
contributions that include the point i located at Xi. Since an 
intersection, and hence a non-zero contribution, is only pos- 
sible if all balls are less than 2r apart, the partial Minkowski 
functionals measure the local morphology in a well-defined 
neighbourhood, namely a ball of radius 2r, around each 
point. 



* Although intersections of four and more poin ts may occur for 
very special point distributions (Platzoder 1995), they are irrele- 
vant since they provide no contribution on average. In practice, 
slight displacements of the order of the numerical resolution can 
be used to resolve such situations without losing accuracy. 



ation of Minkowski functionals for a set of points in a rect- 
angular box with periodic boundary conditions. For detailed 
information, the reader is referred to the source code of the 
program that is available from the authors. Basically, one 
calculates all partial Minkowski functionals a nd su ms them 
up to estimate the density through Equation (A2), by 



1 

-Y^V^A,,,^ 



\D 



(A6) 



Even for a complicated survey geometry, the global 
Minkowski functionals of the Boolean grain model can now 
be written as sums of the partial Minkowski functionals. 
Since only neighbours within 2r around a point contribute 
to its partial Minkowski functionals, we can simply restrict 
the summation to the part of the sample that is further than 
2r from the boundary. Calling this shrunken window D 2r , 
we have 



1 

j— [X)X0 2r (Xi)V M (,<4r;Xi), 



where 



\D 



XD 2r (x) 



1 if x € D 2r 
if x £ D 2r 



(A7) 



(A8) 



is its characteristic function. In the language of spatial statis- 
tics, this quantity is a minus estimator for the volume den- 
sities of the Minkowski functionals. Minus estimators es- 
pecially for the two-point correlation function are already 
known in cosmology; they ha ve been used for example by 
Coleman fc Pietronero (1992) , and thoroughly investigated 
by [Kerscher (1999)| 

One should keep in mind that minus estimators provide 
unbiased estimates only if applied to a sample from a sta- 
tionary point process. Hence we always use volume-limited 
subsamples from catalogues when carrying out our analysis. 



APPENDIX B: DEFINITION OF 
FILAMENTARITY AND PLANARITY 
THROUGH A TOY MODEL 

Consider a random mixture of A^ bodies with different 
Minkowski functionals. The average Minkowski functionals 
of the resulting union are given by (Mecke 1994) 



Vo 



1 



- P v 



Vl = &- pV °pVl 



-pV 



P v 2 - f p>vl 



(Bl) 



V3 



e ~ pVo I pV 3 



-P ViVt + —p v 1 



where V M is the number weighted average over the isolated 
Minkowski functionals of all different bodies. 

Although the Boolean grain model used in our anal- 
ysis places identical balls around all points, we may use 
the formula above for calculating analytically the Minkowski 
functionals of a simple toy model. In fact, we can take into 
account the inhomogeneous distribution of identical balls 
by considering them as randomly distributed objects with 
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Figure Bl. The Minkowski functionals of CfA2N (solid line) 
compared with the best fit of our toy model (dotted line). 



Minkowski functionals characteristic of their respective en- 
vironment. 

Let us assume that points are located in one- 
dimensional filaments, in two-dimensional sheets, or in a 
homogeneous field. If each of these types of structure were 
isolated from the rest, we could calculate the Minkowski 



functionals per single point analytically (Schmalzing 1996) 



For the average point in the filament, sheet and field, re- 
spectively, we haveP] 



. x + e 



vf l) 

Vo(B) 4 x 2 



1-e" 



il 



vf l) 

Vi(B) 



Vf l) _ x 1 - e" 2 * 



(B2) 



V 2 (B) 2 

T7(fll) 
V 3 -2x 

— e 



x I die " v 1 — t 2 arcsint, 
'o 



V 3 (B) 



i We have that ^(x) = -j= J? dtexp(£ 2 ) and 7q(x) is the modi- 



tied Bessel function of order zero. 



—(field) 



0* = o, 



,3) 



(B4) 



V,(B) ' 

where we have normalised the functionals by dividing by the 
values for an isolated ball of the same radius, and x — r/d 
denotes the radius divided by the mean separation of the 
points on the filament or sheet. 

If we further assume that those idealised structures are 
mixed randomly to form the point distribution we wish to 
study, the Minkowski functionals n M are given by Equa- 
tion (Bl), where the V p. are a weighted average over the 



three types of idealised structure. We can therefore calcu- 
late the Vp. numerically, extract the V M , and perform a linear 
fit to obtain the percentage for each type of structure in the 
mixture. By using the quantities 1 — V M /V^(B), which are 
exactly zero for field galaxies, we obtain 



V, 



V„(B) 



T\ 1 



V 



(fii) 



V,(B) 



+ r\ i 



v 



(sheet) 



Vp(B) 



(B5) 



for the toy model. We fix the free parameters T and V by 
minimizing the x 2 OI a ll f° ur Minkowski functional profiles 
at the same time. This procedure yields two numbers for T 
and V\ we interpret them as measures of filamentarity and 
planarity in the point distribution. 

The st andard method o f linear fitting solves the normal 
equations (Press et al. 1987), and also gives an error for our 



measurements, usually of the order of 0.01. As an example, 
Figure Bl compares the Minkowski functionals of the CfA2N 
sample and the best-fitting toy model. Refinements of the 
fitting procedu re, such as a Principal Component Analysis 
(Kendall 1980) are conceivable, but appear inappropriate 
in view of the vast simplifications and resulting systematic 
problems of the toy model. The assumption that galaxies 
arrange in infinitely thin filaments and sheets is of course 
a strong simplification. Nevertheless, this simplified model 
already allows us to put a quantitative measurement on our 
intuitive impression of the local Universe. 
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